Mean field analysis of quantum phase transitions in a periodic optical superlattice 
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In this paper we analyze the various phases exhibited by a system of ultracold bosons in a periodic 
optical superlattice using the mean field decoupling approximation. We investigate for a wide range 
of commensurate and incommensurate densities. We find the gapless superfluid phase, the gapped 
Mott insulator phase, and gapped insulator phases with distinct density wave orders. 

PACS numbers: 03.75.Nt, 05.10.Cc, 05.30.Jp, 73.43Nq 



C/2 



o 



> 

in 

oo 
m 

rn 

O 



X 



I. INTRODUCTION 

Mean-field theory has proved to be a useful tool for the 
analysis of the quantum phase transitions in the lattice 
systems [l|, • The zero-temperature phase diagram of 
the Bose-Hubbard model predicting the superfluid (SF) 
- Mott insulator (MI) transition was first discussed by 
Fisher et al [3]. Jaksch et al ^] suggested the pos- 
sibility of such a transition in an optical lattice loaded 
with ultra cold atoms and it was subsequently observed 
experimentally by Greiner et al in 2002 |;4i]. There are a 
number of reviews on this topic (5|-|7[. Several versions 
of the mean-field theory have been used in the context of 
the ultracold atoms; such as the Bogoliubov approxima- 
tion [8j], the Gutzwiller approach [2| and the mean-field 
decoupling approximation [9|. In the weak interaction 
limit, the Bogoliubov approach is useful. However, it 
is not suitable for the study of the SF-MI phase transi- 
tion since it is valid only for weak interactions. In the 
decoupling approximation, the Bose-Hubbard Hamilto- 
nian is decoupled into single-site Hamiltonians. The re- 
sulting mean-field equation can be solved in two ways; 
either analytically using perturbation theory, or numer- 
ically by diagonalizing the Hamiltonian matrix self con- 
sistently using a convenient basis. The Gutzwiller mean- 
field approach has been used in several papers to study 
the Bose-Hubbard model in quantum lattices [2, [13, llll ■ 
In this paper, we have applied the decoupling approxi- 
mation to a d-dimensional periodic optical superlattice 
with a periodicity of two sites 12]. 

A number of papers on the ultracold atoms in differ- 
ent types of optical superlattices have been published in 
the past few years p^ - [l7j . Experiments on this subject 



have been proposed and carried out in different laborato- 
ries [18-20]. In this context, it is desirable to understand 
the possible phases in different kinds of optical superlat- 
tices. The main purpose of this study is to understand 
the phases in the d-dimensional optical superlattice with 
a periodicity of two sites. For this purpose we use the 
mean-field theory in the the decoupling approximation 
to convert the full Hamiltonian into a sum of single cell 
Hamiltonians. Our findings for ultracold atoms in the 
optical superlattice with a periodicity of 2 sites yields 
gapped insulators accompanied by different crystalline 
orders in addition to the usual Mott insulator and su- 
perfiuid phases. These unusual insulating phases have 
been generically referred to as the superlattice induced 
Mott insulators (SLMIs) in the literature [l3|- Unlike the 
normal Mott insulator phase where strong on-site inter- 
atomic interactions gives rise to the gapped insulator, the 
SLMI phases arise due to the superlattice potential. De- 
pending upon the distribution of bosons within the unit 
cell, there can be various types of SLMI phases. If the 
configuration of the occupancy of bosons within the unit 
cell is such that the alternate sites are occupied by one 
atom and the other being empty, then such an insulator 
is called SLMI-I. The configuration where the alternate 
sites are occupied by two bosons and the other is empty, 
is called SLMI-II. If the configuration is such that alter- 
nate sites are occupied by two bosons, and the other by 
one, then it is called SLMI-HL 
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The rest of the paper is organized in the following man- 
ner. In the next section, we describe the application of 
the mean-field decoupling approximation to an optical 
superlattice with a periodicity of two sites. In section HI 
we present our results. Our conclusions are given in the 
last section; i.e. section IV. 



II. MEAN-FIELD CALCULATIONS FOR THE 
OPTICAL SUPERLATTICE 



which can be written as a sum of single-site Hamiltonian, 
i.e., 7?*^^ = X;,^f^, where 



The system of bosons in a general optical superlattice 
can be best described by the Bose-Hubbard model as 
follows:- 
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In the above equation, (i, j) denotes a pair of nearest 
neighbor sites i and j, t denotes the hopping amplitude 
between adjacent sites, U represents the on-site inter- 
atomic interaction, a\{ai) is the creation (annihilation) 
operator which creates (destroys) an atom at site i and 
hi = a[ai is the number operator, fXi represents the on- 
site chemical potential. For an optical lattice fii = fi for 
all i. However, this is not true for the optical superlat- 
tices and explicit dependence of fj.i on the lattice site i 
depends on the specifics of the superlattices. 

The most important step in obtaining the mean-field 
Hamiltonian is the decoupling of alaj into single site 
operators. For this purpose, we make the following 
approximation:- 

ai = (jji + df, a] = (j)* + a] 

Here 4>i = (a.;), is the mean value and the superfluid 
order parameter, and a is the small fluctuation over the 
mean value. We assume (j)i to be real [l|, hence, 0^ = 0* 
for all i. Substituting the above approximation in the 
kinetic energy term of the Eq. ([1]) , we get 

- t Yl ("1°"^ + '^■'^) = ^*Y1 (^«^^j + ^*s^) 

(id) (iJ) 



-t2_^{dl^j + aj4>i + di(j)j + a](j}i -|- 20^0^) (2) 

We neglect the first term, which is second order in 
fluctuation. The validity of such an approximation can 
be assumed when t is small compared to the the in- 
teraction U and the superlattice potential A. Defining 
0i = -X)(5'^i+<5' ^ being summed over z — 2d nearest 
neighbors with d being the dimension of the optical lat- 
tice, we get the following mean-field Hamiltonian, 



H' 



= -iz y^ [4>i {d\+di) + 4>i 4>i 



-I-— ^ni(ni - I) -^^i^i (3) 



zt 



= -(l)i{a\ + a.i) + (p^cjii + —niifii - 1) - jiiUi. (5) 



We have divided the single-site mean-field Hamiltonian 
by zt to make it and other parameters dimensionless, 
thus U = U / zt, jli = Hi/zt are dimensionless on-site 
interaction and chemical potential respectively. 

For an optical lattice, all the sites are equal, thus 
Hi — IX and (pi — cj) for all i. The Hamiltonian H^^^ 
can then be diagonalized in the following manner. As- 
suming an initial value for the superfluid order parameter 
0, the matrix elements of the mean-field Hamiltonian is 
constructed in the number occupation basis \ n >, where 
71 = 0, 1, 2, • • • , Timax i whcrc Umax IS the maximum num- 
ber of bosons allowed per site whose value depends on 
the on-site interaction U and the chemical potential /x. 
Relatively small values of rimax should suffice for large 
values of U and small values of /i and vice versa. Since 
we have taken four different values of U ranging from 2 
to 15 and the density is always less than 4, we have taken 
iT-rnax = 10 in our Calculation. The Hamiltonian matrix is 
then diagonalized to obtain the lowest eigenstate which 
is used to obtain the new value for 0. Using this value 
for (p, the calculation is repeated till (p is converged. 

We now extend the mean-field for the optical superlat- 
tice which is formed by the superposition of two optical 
lattices with different wavelengths and a relative phase 
shift with respect to each other. This superlattice has a 
periodicity of two sites, thus each unit cell consists of 2'^ 
sites with alternate sites have an energy shift A^. Such 
a system can be described by the Bose-Hubbard model 
([1]) taking into account the relative energy shifts of the 
potential minima such that ^i = /^ — Aj, A^ denotes the 
energy shift and is called the superlattice potential. 

In this optical superlattice, the sites are not equal since 
/ii — /I — Xi is not same for all i. However, the difference 
is restricted only within the unit cell and the whole sys- 
tem is built using this unit cell. The unit cell in the 
superlattice considered here consists of 2'^ sites. Since all 
dimensions are equivalent and the periodicity is two in 
our system, we can work on any one such direction, and 
hence we denote the two neighboring sites by 1 and 2. 
The mean-field Hamiltonian for such a unit cell can be 
written as 



H. 



Substituting di — Ui — (pi in the above equation, we get 



^'^ = -02(al+ai)-0i(4 + a2) + 20102 

+ y[«i(ni - 1) + n2("-2 - 1)] 
-jl[ni + ^2] + Aini + A2n2 



(6) 



H^iP = ^tzY,[Ha\ + a^) - M^] 
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In order to diagonalize the above Hamiltonian, we ex- 
press all the operators including the Hamiltonian in the 
occupation number basis. Then we take initial guess val- 
ues of the superfluid order parameters, 0i and 02. After 
diagonalising the -ff *c^ matrix using the standard Jacobi 
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FIG. 1: (Color online) Variation of average density p as a 
function of the chemical potential /i for U=2, but for different 
values of A starting from 0.5 (red solid curve) to 5.5 (magenta 
double dash dot curve) at the intervals of 1.0. 



FIG. 3: (Color online) Variation of average density of a unit 
cell p as a function of the chemical potential p for U=5, but 
for different values of A starting from 0.2 (red solid curve) to 
7.2 (orange large dashed curve) at intervals of 1.0. 



method, we find the ground state energy, and the ground 
state wave- function. From the relation (/)i = {p-i), we cal- 
culate the superfluid order parameters using the ground 
state wave-function. We then substitute these new val- 
ues of 01 and (^2 in the H^c and iterate the process, until 
the values of (\)\ and (\)2 converge to 10~^. The different 
phases are then analysed based on the values of these 
superfluid densities. 



III. RESULTS 

Taking the superlattice potential for the two distin- 
guished sites within the unit cell Ai = and A2 = A, 
we present our results for a wide range of A, densities p 
and four characteristic values of the on-site interaction 
[/, chosen to cover a substantial part of the phase dia- 
gram. In our analysis we have taken t = 1.0, so all other 
quantities are expressed in units of t. 

First we investiffate the effect of the suDcrlattice do- 




FIG. 2: (Color online) Variation of average superfluid density 
p^ as a function of p for the same set of parameters as in 

Fig.m 



tential on the superfluid phase. In the Figs. [T] and [H we 
plot, respectively, for U — 2 the average density p and 
the superfluid density p* as a function of the chemical 
potential p for different values of A starting from 0.5 to 
5.5 at an interval of 1.0. Here p and p* are, respectively, 
the average density and the superfluid density of a unit 
ceh, i.e., p = (pi + P2)/2, p" — {(j^l + 0i)/2. It is known 
that for U = 2, X — 0, the model ([1]) is always in the su- 
perfluid phase [1| irrespective of the value of the density 

P- 

From the Figs. [T] and [5J we find that density p increases 
with increase in p for finite, but small values of A. The 
superfluid density p^ remains finite for all densities which 
imply that the system continues to be in the superfluid 
phase as in the case of A = 0. However, as A is increased 
further, say for example, A = 4.5, the density develops a 
plateau at p = 1/2 for a range of p values. This is the sig- 
nature of a finite gap in the energy spectrum in this range 
of p values and vanishing of compressibility k = -tA. The 
Fig. [21 suggest vanishing of superfluid density in the same 
range of /i. For all other densities, i.e., p 7^ 1/2, including 
integer densities, the superfluid density remains flnite. 
These features confirm that, for U = 2, model ([T|) is in 
the superfluid phase for all values of p 7^ 1/2 for all val- 
ues of A. However, for p = 1/2 there is a superfluid to an 
insulator phase transition as A is increased. This insu- 
lator phase is different from the standard Mott insulator 
phase arising due to the on-site interaction. Here U = 2 
is small and the Mott insulator phase is not expected. 
The reason for the formation of an insulator phase for 
U = 2 is due to the superlattice potential, and to dis- 
tinguish this insulator from the Mott insulator phase, we 
call it as superlattice induced Mott insulator (SLMI) [13| 
as mentioned earlier. In order to understand the SLMI 
phase, the distribution of bosons within the unit cell is 
tabulated in the Table ID As we discussed in the previ- 
ous section, the unit cell consists of 2^* sites and each cell 
has two distinct sites, which we refer to as 1 and 2. The 
values of site densities pi , p2 and superfluid densities pf , 



P2 are listed in the table for different values of A. For 
A < 3.7, the on-site superfluid densities p\ and P2 remain 
finite for all densities. However, for A > 3.7 and density 
p — 1/2, we find p\ — P2 — ^ and pi ~ 1, p2 = 0. This 
implies, within the unit cell, one site is occupied and the 
other site being empty. Since this unit cell repeats to 
cover the entire lattice, it has every alternate site occu- 
pied and the other is empty like a charge density wave 
(CDW) phase which normally arises due to the nearest 
neighbor interaction. However, it should be noted here 
in this work that the CDW like density distribution is 
due to the superlattice potential and there is no nearest 
neighbor interaction involved. Since the distribution of 
bosons follow a pattern [10 10 10 • • • ] in all the 
d directions of the lattice, we call this phase as SLMI-I 
to distinguish it from SLMI-H discussed below. The Ta- 
ble U] also confirms that there is no insulating phase for 
densities p = 1 and 2. 

Results for U = h are similar to those of [/ = 2. In 
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FIG. 4: (Color online) Variation of average superfluid density 
of a unit cell as a function of ^ for the same set of parameters 
as in Fig. [S] 





TABLE l:U = 2.0 




A 


P 


M 


Pi 


Pi 


m 


n2 


0.5 


0.620 


0.00 


0.58 


0.45 


0.73 


0.51 


0.5 


1.022 


0.70 


0.87 


0.75 


1.12 


0.92 


0.5 


2.017 


2.60 


1.72 


1.57 


2.13 


1.91 


1.5 


0.490 


0.20 


0.52 


0.21 


0.77 
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1.020 
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0.98 


0.62 


1.32 


0.71 
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2.020 
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0.47 


0.39 
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0.91 


0.08 
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0.49 
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FIG. 5: (Color online) Variation of average density of a unit 
cell p as a function of the chemical potential p for U=10, but 
for different values of A, varying from 0.2 (red solid curve) to 
14.2 (or " " ■ ' ' 




FIG. 6: (Color online) Variation of average superfluid density 
of a unit cell as a function of p for the same set of parameters 
as in Fig. (5] 



Figs. [3] and |4] we plot, respectively, density p and the 
superfiuid density p^ as a function of p for different values 
of A. The system is in the superfiuid phase a.t p — 1/2 
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A 
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Pl 
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ni 
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0.523 


0.0 
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0.35 


0.57 


0.48 
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0.270 
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FIG. 7: (Color online) Variation of average density of a unit 
cell p as a function of the chemical potential /i for U=15, but 
for different values of A, varying from 0.2 (red solid curve) to 
18.2 (violet large dot dashed curve) at intervals of 3.0 
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FIG. 8: (Color online) Variation of average superfiuid density 
of a unit cell as a function of [i for the same set of parameters 
as in Fig. [T] 



and p — \ initially for low values of A (< 2.6). For 
A > 2.6, a plateau appears in the p versus p, plot for 
p = 1/2 suggesting a gap in the energy spectrum. The 
superfluid density vanishes in this region. This plateau 
at p = 1/2 widens as A increases. However, the system 
remains in the SF phase at p = 1 for all the values of 
A considered. In Table |TT1 we tabulate the values of site 
densities and superfluid densities within the cell and we 
conclude that the transition from the SF to the SLMI-I 
phase is at A = 2.6, when the superfluid density vanishes, 
and the occupancy configuration is of the form [10 10 
• • • ] . On the other hand, at other values of p and for all 
values of A, the superfiuid densities, p\ and p^^ remain 
finite. 

Results for [/ = 10 are different from that oiU — h. 
This difference is mainly due to the fact that model ^ 
has SF to MI transitions for integer densities. For p = 1, 
the critical Uc ~ 5.8 for the SF-MI transition [ij and 
this implies that for U = 10, model ([T]) is in the Mott 
insulator phase for p = 1. In the Figs. [5] and [6] we plot. 



TABLE III: V = 10.0 
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Pi 


P2 


ni 


n2 


0.2 


0.48 


0.00 


0.30 


0.29 


0.53 


0.44 


0.2 


1.00 


2.00 


0.00 


0.00 


1.00 


1.00 


0.2 


2.00 


14.50 


0.00 


0.00 


2.00 


2.00 


2.2 


0.50 


1.20 


0.00 


0.00 


1.00 


0.00 
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1.00 
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0.10 


0.00 


0.00 


1.00 
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0.00 
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0.00 
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1.00 



respectively, density p and the superfiuid density p^ as a 
function of p for different values of A. From these figures 
the following conclusions are drawn. For small values of 
A, the plateau in the p versus p plot exist only for p = 1 
and the p* vanishes in the same range of p confirming 
the expected SF to MI transition for p = 1. The system 
remains in the superfiuid phase for all other densities. 
However, as we increase A, the plateau region at p = 1 
, i.e., the MI region, shrinks first, completely disappears 
for some values of A, and re- appears again for higher 
values of A. 

A plateau develops at p = 1/2 for A > 2.3 and at 
p = 3/2 for A > 5.3. From the tabulated values of pi 
and p2 in the Table iHll the insulator phase at p = 1/2 
is the same as SLMI-I. The insulator phase at p = 3/2 
has a density distribution [2 1 2 1 ■ • • ] which we call 
SLMI-IH. The insulator phase at p = 1 for higher values 
of A has the occupation at alternate lattice sites [2020 
• • • ] which we refer to as the SLMI-II phase. 

Thus for U — 10.0, for A < 2.3 the system exhibits a 
Mott insulator phase for p = 1 and SF phase else where. 
For 2.3 < A < 5.3, the system has two insulating phases: 
SLMI-I for p = 1/2 and MI phase for p = 1. The system 
is in the superfiuid phase for the rest of the densities. 
For A > 5.3 the system shows SLMLI for p = 1/2, MI 
phase for p = 1, SLMLIII for p = 3/2 and SF for other 
densities. For p = 1 the MI phase is lost for A > 6.5 
and re-appears as SLMI-II for A > 13.1. The results 
for U = 10.0 are in qualitative agreement with those 
obtained using DMRG [13 ■ 

The system for U = 15.0 behaves similarly as f/ = 10.0. 
At p = 1/2, the system starts off in the gapless SF phase 
for low values of A(= 0.2), as evident from the Figs. [7] 
and [5] and Table HVl But at p = 1.0 and 2.0, the system 
is in the MI phase at this value of A. As A is increased 
to a value greater than 2.2, a gap appears at p = 1/2, 
marking the transition from the SF to the SLMI-I phase, 
as seen in the Table IIVI where we have vanishing su- 
perfiuid densities, and also a occupancy configuration 
of [1 1 •■•]. Also, at p = 3/2, another gap ap- 



pears at A = 4.8, implying the transition from the SF 
phase to the gapped SLMI-III phase with a configura- 
tion [2 f 2 f • • •]. As A becomes greater than 11.7, the 
system at p = 1.0 undergoes a phase transition from the 
MI to the SF phase, shown by the non-zero values of the 
superfluid density. As A becomes larger than 18.1, the 
gap reappears once again, showing that the system has 
entered into the gapped SLMI-II phase with configura- 
tion [2020 •••]. 



IV. CONCLUSIONS 

We have analyzed the various phases exhibited by a 
system of bosons in an optical superlattice with a unit 
cell consisting of two distinct lattice sites using the mean- 
field decoupling approximation, for various values of the 
superlattice potential. A, corresponding to four values of 
the on-site interaction U. For U = 2.0, we find that the 
system resides in the SF phase for all densities for small 
values of A. At p = 1/2, there is a transition from the SF 
to the SLMI-I phase at A ~ 3.7, but for other densities, it 
remains in the gapless SF phase. For U = 5.0, the system 
undergoes a SF - SLMI-I phase transition for p = 1/2 at 
A = 2.6, but remains in the SF phase for other densities 
and A. For U = 10.0, the system undergoes a SF - SLMI- 
I phase transition at A = 2.3 for p = 1/2,. However, for 
p = 1, the system starts in the MI phase, as the value of 
U is large, and as A is increased, the gap in the MI phase 
shrinks, and eventually goes to zero, marking the MI-SF 
phase transition at A = 6.5. The system stays in the gap- 
less SF phase for 6.5 < A < 13.1. As A is increased further 
the system undergoes a phase transition from SF - SLMI- 
II at A = 13.1. For p ~ 3/2, we see a phase transition 
from SF to SLMI-III phase at A = 5.3. Similar behavior 
is observed for U = 15.0, with the system for p = 1/2 
undergoing a phase transition at A — 2.2. For p = 1, the 
system makes a transition from the MI to the SLMI-I 
phase at A = 11.7, and then from the SF to the SLMI-II 



phase at A = 18.1. Also a phase transition from the SF 
to the SLMI-III phase at A = 4.8. It should be possible 
to extend this calculation to superlattices with different 
periodicity. The charge density wave order in the SLMI 



TABLE IV; U=15.0 



A 


P 


M 


Pt 


Pi 


111 


n2 


0.2 


0.52 


0.10 


0.28 


0.28 


0.57 


0.47 


0.2 


1.00 


1.30 


0.00 


0.00 


1.00 


1.00 


0.2 


2.00 


17.80 


0.00 


0.00 


2.00 


2.00 


6.2 


0.50 


0.18 


0.00 


0.00 


1.00 


0.00 


6.2 


1.00 


6.66 


0.00 


0.00 


1.00 


1.00 


6.2 


1.50 


16.11 


0.00 


0.00 


2.00 


1.00 


15.2 


0.50 


0.10 


0.00 


0.00 


1.00 


0.00 


15.2 


1.02 


15.10 


0.60 


0.30 


1.55 


0.49 


18.2 


0.50 


0.09 


0.00 


0.00 


1.00 


0.00 


18.2 


1.0 


16.00 


0.00 


0.00 


2.00 


0.00 


18.2 


1.50 


19.08 


0.00 


0.00 


2.00 


1.00 



phase will depend on the number of distinct sites within 
the unit cell. The mean field approach is exact in the in- 
finite dimension and the error, because of neglecting the 
fluctuations, become severe in low dimensions 2l|. How- 
ever, it proves to be an excellent tool for the qualitative 
analysis (e.g. phase diagram), which is our focus in this 
paper. Since the parameters of the Hamiltonian can be 
varied to a large range of values by tuning the strength 
of the optical potentials, we hope our detailed study of 
model ^ will stimulate experimental studies that could 
lead to the observation of the superlattice induced Mott 
insulators. 
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